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The eigen-frequencies of the axial w-modes of oscillating neutron stars are studied using the 
continued fraction method with an Equation of State (EOS) partially constrained by the recent 
terrestrial nuclear laboratory data. It is shown that the density dependence of the nuclear symmetry 
energy Eaym{p) affects significantly both the frequencies and the damping times of these modes. 
Besides confirming the previously found universal behavior of the mass-scaled eigen-frequencies as 
functions of the compactness of neutron stars, we explored several alternative universal scaling 
functions. Moreover, the u;//-mode is found to exist only for neutron stars having a compactness of 
M/R > 0.1078 independent of the EOS used. 
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I. INTRODUCTION 

Presently, among the greatest challenges in modern 
physics are determining the equation of state (EOS) of 
dense neutron-rich matter and detecting gravitational 
waves - tiny ripples in space-time predicted by the theory 
of general relativity. Physics of neutron stars provides a 
unified framework for studying both phenomena: on the 
one hand neutron stars are natural sites of super-dense 
matter, while on the other hand they are one of the major 
candidates for emitting gravitational waves in the band- 
width targeted by the ground based laser interferometric 
detectors LIGO VIRGO (e.g. 0) and GEO Q. Be- 
sides the gravitational waves from elliptically deformed 
pulsars, an important mechanism of gravitational radia- 
tion from a neutron star is its non-radial oscillations. The 
latter may provide potentially a unique probe to the EOS 
of super-dense neutron-rich matter because the resulting 
gravitational waves carry important information about 
the neutron star structure [1,1^, Q . The non-radial neu- 
tron star oscillations could be triggered by various mech- 
anisms such as gravitational collapse, a pulsar "glitch" 
- a sudden disrupt of the otherwise very regular pulsar 
period, gravitational spin-orbit coupling due to a com- 
panion in a close binary orbit, or a transition to novel 
phases of matter in the inner core. 

In the framework of general relativity, gravitational ra- 
diation damps out the neutron star oscillations. The 
frequency of the non-radial oscillations thus becomes 
"quasi-normal" (complex) with a real part representing 
the actual frequency of the oscillation and an imagi- 
nary part representing the energy losses. The eigen- 
frequencies of the quasi-normal modes can be found by 
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solving the equations describing the non-radial pertur- 
bations of a static neutron star in general relativity 
0, @|- These equations are derived by expanding the 
perturbed tensors in spherical harmonics. The pertur- 
bation equations split into two sets: one for the axial 
mode for which the spherical harmonics transform un- 
der parity (i.e., under the transformation — > tt — and 
(j) ^ TT -\- (j) ) a.s (—1)'+^; and second for the polar mode 
for which the spherical harmonics transform under parity 
as (—1)', respectively. With the appropriate boundary 
conditions at the neutron star center and surface, the so- 
lution of the perturbation equations yields the complex 
eigen-frequencies. They can be classified as follows: (1) f- 
modes, associated with the global oscillation of the fluid; 
(2) g-modes, associated with the fluid buoyancy; and (3) 
p-modes, associated with the pressure gradient. These 
modes exist in both Newtonian gravitational theory and 
general relativity. There is still another type of mode that 
only exists in general relativity - the so called w-mode 
associated with space-time and for which the motion of 
the fluid is negligible [1, 0, [1| • The w-mode is very im- 
portant for astrophysical applications since it is related 
to the space-time curvature and exists for all relativistic 
stars, including black holes. The standard axial w-mode 
is categorized as wj. Additionally, there exists an inter- 
esting family of axial w-modes [9|, categorized as wij. 
According to the work of Chandrasekhar & Ferrari 0, d] , 
the axial w-mode is described by a unique second-order 
differential equation. Its simple form makes it possible 
to conduct detailed numerical studies of its spectrum. A 
major characteristic of the axial w-mode is its high fre- 
quency accompanied by very rapid damping. In fact, the 
w-mode frequencies are outside the peak of the detec- 
tor sensitivities of the currently operating and planned 
gravitational wave detectors. Nevertheless, it is theoret- 
ically interesting to investigate the w-mode of neutron 
star oscillations as it is an important mechanism for pro- 
ducing gravitational waves 
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study will help stimulate discussions on detecting high- 
frequency gravitational waves with new generations of 
detectors in the future. 

The gravitational wave frequency of the axial w-modc 
depends on the structure and properties of neutron stars 
, which arc determined by the EOS of neutron-rich stel- 
lar matter. However, at present time the EOS of mat- 
ter under extreme conditions (densities, pressures and 
isospin asymmetries) is still rather uncertain and theo- 
retically controversial. One of the main sources of uncer- 
tainties in the EOS of neutron-rich matter is the poorly 
known density dependence of the nuclear symmetry en- 
ergy, Esym{p) The EOS of neutron-rich nuclear 
matter plays an important role not only in many astro- 
physical objects/processes [31 but also in heavy-ion colli- 
sions especially those induced by neutron-rich radioactive 
beams in terrestrial laboratories While heavy -ion 
collisions are not expected to create the same matter and 
conditions as in neutron stars, the same elementary nu- 
clear interactions are at work in the two cases [15| . Thus, 
it is important to examine ramifications of conclusions 
regrading the EOS extracted from one field in the other 
one. On one hand, extremely impressive progress has 
been made in astrophysical observations relevant for con- 
straining the EOS of nuclear matter. To our best knowl- 
edge, however, mainly because of the low precision associ- 
ated with the current measurements of neutron star radii, 
a non-controversial conclusion on the EOS and the den- 
sity dependence of the nuclear symmetry energy has yet 
to come from analyzing astrophysical observations. Nev- 
ertheless, it is very interesting to note that since the pio- 
neering work of Lindblom [16[ , a lot of efforts have been 
devoted to extracting the underlying EOS by using the 
technique of inverting the TOV (Tolman-Oppcnheimer- 
Volkov) equation using the masses (M) and radii (i?) 
simultaneously measured accurately for several neutron 
stars, see, e.g., refs. [13, 113 j for the latest reports. While 
the principle of this technique is well demonstrated and 
promising, the lack of simultaneously measured (Af , R) 
data has so far hindered the fruitful applications of this 
technique. In particular, while the masses in some bi- 
nary systems, such as the double neutron star binaries, 
are well measured, their radii arc unfortunately not pre- 
cisely known [l3, EB] • Similarly, under the premise that 
the frequencies and damping rates of several w-modes 
of a single neutron star are measured, Tsui and Leung 
[l9| investigated an inversion scheme to determine the 
mass, radius and density profile of the neutron star. The 
underlying EOS can then also be inferred. While the ap- 
proach is convincingly promising, it can only be applied 
after the gravitational wave astronomy becomes a real- 
ity. On the other hand, heavy-ion reactions especially 
those induced by radioactive beams provide an alterna- 
tive means to constrain the EOS of neutron-rich nuclear 
matter, see, e.g., refs.[il [13, (H HI HI, for reviews. 
Over approximately the last 40 years, the nuclear physics 
community has made steady progress in constraining the 
EOS of symmetric nuclear matter up to several times 



the normal nuclear matter density. However, it was only 
during roughly the last 10 years the nuclear physics com- 
munity has made significant progress in constraining the 
density dependence of the symmetry energy of neutron- 
rich nuclear matter thanks mainly to the rapid technical 
advantages in accelerating radioactive beams [iJl- It is 
particularly exciting to see that some significant results 
in constraining the density dependence of the symm etry 
energy have beeb reported very recently [13, [H, [1^, [27| . 
Thus, it is interesting to investigate how the EOS con- 
strained by the available data from heavy-ion reactions 
may help limit the ranges of frequencies and damping 
times of the w-mode and the strain of other forms of 
gravitational waves [2^. In this work we apply an EOS 
with its symmetric part constrained up to about 5 times 
the normal nuclear matter density and its symmetry en- 
ergy constrained in the subsaturation density region by 
the available heavy-ion reaction data. We have two major 
objectives: (1) to investigate to what extent the nuclear 
symmetry energy could affect the eigen-frequency of ax- 
ial w-mode; and (2) to examine the relationship between 
the axial w-mode eigen-frequency and the neutron star 
gravitational field. This study is important because the 
axial w-mode is believed to be a pure space-time mode, 
its eigen-frequency is expected to have a more direct re- 
lationship with the properties of gravitational field, e.g. 
with the energy of gravitation. Compared to previous 
studies on the w-mode in the literature, the major issues 
we address here arc scientifically important and studied 
from a different direction. In light of the great impor- 
tance of and the extreme difficulties encountered in de- 
tecting gravitational waves, our results presented here are 
scientifically useful. 

This paper is organized as follows: after the introduc- 
tion, in Section II we recall the formalism and the contin- 
ued fraction method for computing the eigen-frequency 
of the axial w-mode (the set of perturbation equations 
and their numerical solution); we then briefiy outline the 
EOSs applied in our studies in Section III; in Section IV 
we present our numerical results; and in Section V we 
conclude with a short summary. 



II. FORMALISM FOR CALCULATING THE 
AXIAL W-MODE FREQUENCY 

In what follows we summarize the formalism for calcu- 
lating the eigen-frequency of the axial w-mode and the 
continued fraction method. Unless noted otherwise, we 
use geometrical unit (G — c = 1). According to Chan- 
drasekhar & Ferrari [3, [^ , the axial perturbation equa- 
tions for a static neutron star can be simplified by intro- 
ducing a function z(r), constructed from the radial part 
of the perturbed axial metric components. It satisfies the 
following differential equation 



drl 



+ [c^2 _ v{r)]z = 0, 



(1) 
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where u>{= ujo+iuii) is the complex eigen-frequency of the 
axial w-mode. Inside the star, the tortoise coordinate 
and the potential function V are defined by 



where y = l — a/r,R<a< 2R {R is the stellar radius) 
and the coefficients a„ satisfy the following four-term re- 
currence relation 



dr (or 



d 



dr 



and 



V = -^[1(1 + l)r + 4nr^(p - p) - 6m] 



(2) 



(3) 



with I the spherical harmonics index (used in describing 
the perturbed metric, here only the case I = 2 is consid- 
ered), p and p the density and pressure, and m the mass 
inside radius r, respectively. The functions e"^ and are 
given by the line clement for a static neutron star as 



ds^ 



e^^dt^ + e^^dr^ + r\d0'^ + sin^edcf)^). (4) 



Outside the neutron star. Eqs. [2] and [3] reduce to 

= 2Mln(r - 2M) {or = r - 2M d ^ 
dr^ r dr 



and 



V 



2M 



+ l)r ~ 6M]. 



(6) 



where M is the total gravitational mass of neutron 
star. The solutions to this problem are subject to a 
set of boundary conditions (BC) constructed by Chan- 
drasekhar & Ferrari 0] - regular BC at the neutron star 
center, continuous BC at the surface and behaving as a 
purely outgoing wave at infinity. 

As discussed in references [1, [13] , the w- modes have 
larger imaginary parts and this introduces considerable 
difficulties to dealing with the boundary conditions at 
infinity. For a larger radius, the term representing the 
incoming wave, z*" ^ e"""*"*, vanishes and therefore the 
problem no longer satisfies the necessary BC. A practical 
way to overcome these difficulties is to utilize the contin- 
ued fraction method 0, [sot , which is applicable to both 
the axial and polar w- modes. This method provides both 
the wj- and w/7-modes. 

Here we give a brief introduction to the continued frac- 
tion method following closely Benhar et al. [§| . It is con- 
venient to use the dimensionless geometrical units, i.e. 
2M = c = G = 1. First, one represents the solution 
outside the neutron star as 



+ /Snftn + 7nan-l + l5nan-2 = 0; (n > 2) (8) 

with 



Q!„ = (1 )n{n + 1), 

a 



, 3n 
Pn = -2{iuja + n- —)n, 
Za 



(9) 



(10) 



7n = (l--)"(^-l) + --^(^ + l), (11) 

a a 



Sn = -{n-3){n + l). 



(12) 



The coefficients oq and ai are determined by the conti- 
nuity of z(r) and z(r)^r at r = a, that is 



a-o 



X(a)' 



(13) 



ai 



^ a + T^ a • 14 

X(a) a - 1 



Since Eq. [5] is a four-term recurrence relation, to deter- 
mine a„ uniquely one must know three initial terms, i.e. 
Egs. [T3l and [T4l are not sufficient alone. Leaver [sij has 
shown that Eq. [8] can be further reduced to a three-term 
recurrence relation as 



+ /Jnftn + 7ria„_l = 0, {n > 2) , (15) 



for n = 



for n — 1 



for n > 2 



do = -1, $0 



01 

ao ' 



di = ai, $1 = /3i, 71 = 71, 



(16) 



(17) 



an = an, Pn= Pn T , 7n = In T ■ (18) 



7n-l 



7n-l 



z(r) — (r — ly 



E 



any ' 



According to Leaver [31j and Wall [34I, the function z{r) 
(7) describes a purely outgoing wave at infinity only if 
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/n(w) = /3„ 



-i7« 



-27r. 



/30 



Q!K7n+l 

an+l7n 



(n = 0,l,2,...). 



(19) 



For n = 0, the above expression becomes 

/oM^/3o- =0. (20) 

Eq. [19] is the key expression for computing the eigen- 
frequencies of the axial w-modes. Obviously fa{i-u) = 
is the simplest case, which is mainly employed in our 
calculation. In order to rule out the pseudo roots, the 
cases n = 1 and 2 are also examined. 

The numerical procedure goes as follows: (1) choose a 
suitable range for ujq and w^; (2) integrate Eq. [1] starting 
at the center to r = a; (3) compute /o(<^) through Eq. 
1201 and plot the curve loq against uji (where the real part 
and the imaginary part of function fo{w) are zero); (4) 
determine the crossing point of the curves - the crossing 
point specifies the eigen-frequcncy of the axial w-mode. 
The same procedure is also repeated for the cases with 
n = 1 and n = 2. 



III. THE EQUATION OF STATE OF 
NEUTRON-RICH NUCLEAR MATTER 
PARTIALLY CONSTRAINED BY RECENT 
TERRESTRIAL NUCLEAR LABORATORY DATA 

The necessary input in solving the Eq. [l]is the EOS. 
While at high densities reached in the core of neutron 
stars, other particles and new phases of matter may ex- 
ist, in this work we consider the simplest model of neu- 
tron stars consisting of neutrons, protons and electrons 
{npe) in beta-equilibrium. For isospin asymmetric nu- 
clear matter, various theoretical studies have shown that 
the energy per nucleon can be well approximated by 

E{p, 6) = E{p, S = 0) + E,y„,{p)S^ + 0{S^) (21) 

in terms of the baryon density p ~ pn + Pp, the isospin 
asymmetry d = {pn - Pp)/{pn + Pp), the energy per nu- 
cleon in symmetric nuclear matter E{p, 6 = 0), and the 
bulk nuclear symmetry energy Egym (p) [H, [H, [H, . 
The corresponding pressure of the npe matter is given by 

= p2 [E'{p,6^Q) + Ei^,MS''] 

+ ^6{l-S)pE,y„,{p), (22) 

where pe = 5(1 - S)p and pe = - Mp = 4(5-Bsy,n(p) 
are, respectively, the density and chemical potential of 



electrons. It is obvious that the pressure is dominated 
by the symmetry energy term in the npe matter near the 
saturation density where the E'{p, (5 = 0) term vanishes. 
The value of the isospin asymmetry 6 at (3 equilibrium 
is determined by the chemical equilibrium and charge 
neutrality conditions, i.e., 6 ~ 1 ~ 2xp with 



Xp « 0.048 [E,y^{p)/E,y^{po)f (p/po)(l ' 2xp)\ (23) 

It is seen that the proton fraction Xp is uniquely de- 
termined by the density dependence of the nuclear sym- 
metry energy Esymip)- It is well-known that while the 
maximum mass of neutron stars is determined mostly by 
the stiffness of the symmetric nuclear EOS E{p,5 = 0), 
the radii of neutron stars are mostly determined by the 
slope of the Egymip) of the most uncertain 

part of the EOS of neutron-rich nuclear matter is the 
density dependence of the Esymip) especially at supra- 
saturation densities [l3|- In the previous studies of 
the eigen- frequencies of the w-mode using various EOSs 
@, Hi IHj [13 1 the focus was on exploring effects of the 
stiffness of the EOS and understanding the scaling be- 
havior of the eigen- frequencies. It is still not clear what 
arc the effects of the Egymip)- In this study we use EOSs 
with the same incompressibility K at normal density but 
different Esym{p)- Our study is thus complementary to 
the existing studies [1, d, [H, [131 • Moreover, the EOS 
of neutron-rich nuclear matter also plays an important 
role in heavy-ion collisions especially those induced by 
neutron-rich radioactive beams in terrestrial laborato- 
ries. Significant progress has been made recently in con- 
straining the EOS of neutron-rich nuclear matter using 
heavy-ion experiments, see, e.g., refs. [3, [13, [13 • In 
particular, experimental data on collective flow and kaon 
production in relativistic heavy-ion collisions have put a 
strong constraint on the EOS of symmetric nuclear mat- 
ter at densities up to about 5 times the normal nuclear 
matter density |22l |. On the other hand, the Esym{p) 
in the sub-saturation density region has also been con- 
strained recently by analyzing data on isospin diffusion 
[H, [13, [13, [m and isoscaling [13, [13] in heavy-ion re- 
actions at intermediate energies. In this work, we use 
the MDI (Momentum Dependent Interaction) EOSs ob- 
tained from Hartree-Fock calculations using a modified 
Gogny interaction [i^]. The parameters of the interac- 
tions are adjusted such that the resulting EOSs satisfy 
all of the above constraints from heavy-ion reactions. 

For the modified Gogny MDI interaction, the baryon 
potential energy density can be expressed as [i^ 
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Au{x)pnPp _|_ Ai{x) 

Po 2po 



(pI + pI) + 



B 



a + l pZ 



.3 I fr{r,p)fT'{r,f) 



(24) 
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FIG. 1: (Color online) The symmetry energy (a) and pressure 
(b) versus the baryon number density (baryon number density 
p is in units of saturation nuclear density po)- 



In the above equation, the isospin t = 1/2 (—1/2) is 
for neutrons (protons). The cocfhcicnts Au{x) and Ai{x) 
are 



Auix) 



-95.98- 



2B 



1' 



Aiix) 



-120.57 + a;- 



2B 



<J+l 
(25) 

The values of the parameters arc a = 4/3, B = 106.35 
MeV, Cr,T = -11.70 MeV, C^.-r = -103.40 MeV and 
A = which is the Fermi momentum of nuclear matter 
at pq. The parameter x was introduced to mimic various 
Esyrnip) predicted by different microscopic many-body 
theories. By adjusting the x parameter, the Esym{p) 
is varied without changing any property of symmetric 
nuclear matter and the symmetry energy at saturation 
density as the x-dependent A„(x) and Ai{x) arc auto- 
matically adjusted accordingly. We note especially that 
the symmetry energy at normal density Esymipo) is fixed 
at 31.6 MeV and the incompressibility K of symmetric 
nuclear matter at normal density is fixed at 211 MeV 
consistent with known experimental constraints [T3 | . We 
stress that both the Egy^n^po) and K are independent of 
the X parameter. It was shown that [4l|,|4^ only values of 
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FIG. 2: (Color online) Mass-radius relation for static neutron 
stars. 



X in the range between —1 (MDIx-1) and (MDIxO) are 
consistent with the isospin-diffusion and isoscaling data 
as well as the available measurements of the neutron- 
skin thickness of ^°*Pb. These data only constrain the 
Esyrnip) sub-saturation densities. At supra-saturation 
densities, some experimental indications on the trend of 
the Esyrnip) -ti^s just start to appear For this ex- 
ploratory study we assume that the two limiting func- 
tions of Esyrnip) with a; = and x = —1 experimen- 
tally constrained only at sub-saturation densities can be 
extrapolated to supra-saturation densities according to 
the MDI predictions. Hopefully, future theoretical and 
experimental studies will allow us to better constrain 
also the high-density symmetry energy and thus relax 
the above assumption. Moreover, as we indicated ear- 
lier, extrapolating hadronic EOS to high densities with- 
out considering the possible hadron-QGP (Quark-Gluon- 
Plasma) phase transition is always a problem for neu- 
tron star models. Given our limited knowledge about 
the hadron-QGP phase transition and other properties 
of super-dense matter, for our purpose of studying the 
imprints of nuclear symmetry energy on the w-mode, we 
stick to the simplest model of neutron stars made of the 
npe matter. 

For comparisons, we also employ the widely used EOS 
by Akmal, Pandharipande and Ravenhall (APR) [i^ . 
Shown in Fig. 1(a) is the nuclear symmetry energy con- 
sidered in this work. It is interesting to note that the 
APR prediction on the Esymip) at sub-saturation den- 
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TABLE I: Saturation properties of the symmetric nuclear 
matters EOS. The first column identifies the EOS. The re- 
maining columns presenting the following quantities at nu- 
clear saturation density: saturation baryon density; energy 
per particle; incompressibility; nucleon effective mass; sym- 
metry energy. 

EOS poifm,-'') EsjMeV) K{MeV) m{MeV/d') E,y^{MeV) 



MDI 0.160 
APR 0.160 



-16.08 
-16.00 



211.00 
266.00 



629.08 
657.25 



31.62 
32.60 



sities lies right between the MDIxO and MDIx-1 predic- 
tions. At supra-saturation densities it follows closely the 
MDI Esym{p) with X At nuclear densities below ap- 
proximately 0.07 fm~'^, we supplement the EOSs by those 
of Refs. [43, [3 which are more suitable for the neutron 
star crust. The saturation properties of symmetric nu- 
clear matter for the EOSs used here are summarized in 
Table I. We emphasize that the APR predicts a signifi- 
cantly stiffer incompressibility of 266 MeV compared to 
the MDI EOS which has an incompressibihty of 211 MeV 
independent of the symmetry energy parameter x. 

Shown in Fig. 1(b) is the pressure of the npe matter 
as a function of the baryon density. As one expects, the 
stiffer symmetry energy with a; = — 1 leads to a higher 
pressure compared to a; = 0. As a result, the radii of neu- 
tron stars are consistently larger with a; = — 1 as shown 
in Fig. 2. We notice that the MDIxO and MDIx-1 EOSs 
lead to about the same maximum mass for neutron stars. 
This is because they have the same incompressibility K 
of 211 MeV for symmetric nuclear matter at normal den- 
sity. The APR EOS leads to a significantly larger maxi- 
mum mass because of its stiffer incompressibility but sim- 
ilar radii as the MDIxO prediction because of the similar 
symmetry energy density functionals near the saturation 
density. We notice here that the range of radii between 
the MDIxO and MDIx-1 predictions is much smaller than 
that spanned by the various EOSs considered in ref. p^ . 
In the latter the EOSs are different not only in their sym- 
metry energy functionals but also the incompressibilities 
for symmetric nuclear matter. 



IV. RESULTS AND DISCUSSION 

In this section we present our numerical results for the 
first axial w-mode (z/;/-modc), the w/z-mode, the second 
w-mode (u)/2-modc), and some of the third axial w-modes 
( w/a-modcs). The calculation is performed applying the 
continued fraction method together with the EOSs we 
discussed in the previous Section. 

We first examine effects of the symmetry energy on the 
frequencies and damping times of the w-mode. Shown in 
Figs. 3 and 4 are the frequency and damping time of 
the wi-mode (a) and wn-vaode (b) respectively, as func- 
tions of the neutron star mass. These figures establish 
the relationship between the expected frequencies of the 
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FIG. 3: (Color online) Frequency of wi-moAe (a) and wn- 
mode (b) as a function of the neutron star mass M. 



axial w-modes, for a given EOS, and the stellar mass. It 
is interesting to notice, in Fig. 3, that there is a clear 
difference between the frequencies calculated with the 
MDIxO EOS and those with the MDIx-1 EOS. Since the 
major difference between these two cases is the density 
dependence of the nuclear symmetry energy, it is obvious 
that the symmetry energy has a clear imprint on the fre- 
quencies. Therefore, the symmetry energy constrained 
by nuclear reactions in terrestrial laboratories can help 
determine the expected frequency of the axial w-mode 
and its damping time. 

In Figures 5, 7 and 9 we show the real (a) and imag- 
inary (b) parts of the eigen- frequency of w/-, ui/2- and 
w/z-modes scaled by the mass M as a function of the neu- 
tron star compactness M/R, respectively. These results 
suggest that the scaled eigen-frequency exhibits a univer- 
sal behavior independent of the EOS used as a function of 
the compactness parameter. Similar universal behaviors 
were first found for the polar w-mode by Andersson et 
al. and later for the axial w-mode by Benhar et al. [4§| 
and Tsui et al. [s^, [s^l . As it was pointed out earlier in 
the above references, provided that the masses and radii 
of neutron stars are known, the universal scaling behav- 
iors allows an accurate determination of the w-mode fre- 
quency and damping time of gravitational waves. This is 
very important for guiding the gravitational wave search. 
On the other hand, if both the frequency and damping 
time for a given neutron star are known this could pro- 
vide information on the neutron star mass and radius. In 
Figs. 5 and 7 we also display the curves (Fit) best fitting 
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FIG. 6: (Color online) Eigen-frequency, lo, of the w/-mode 
scaled by the gravitational energy \Eg\ versus the metric func- 
tion e^^ at the stellar surface. 
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the results obtained by Tsui et al. using eight different 
EOSs [36|. It is seen that our numerical results are in 
good agreement with theirs. 

Realizing that the axial w-mode is a space-time mode 
0, it is interesting to study whether the w-mode eigen- 
frequencies may be scaled by the total gravitational en- 
ergy and how the scaled frequencies behave as functions 
of the metric function e^^^^) ^ 1/(1 - 2Af/i?) (at the 
stellar surface). This approach relates again the eigen- 
frequency to the compactness parameter M/R of neutron 
stars, but additionally it has the advantage of taking into 
account more completely the space-time properties. The 
gravitational energy is calculated from [50| 



Eg = r A7:r\l - [1 - '^^]-y^}pdT 
Jo r 



(26) 



Shown in Figs. 6, 8 and 10 are the real (a) and imagi- 
nary (b) parts of the eigen- frequencies of the wi-, wi2- 
and tu/z-modes scaled by the absolute value of the grav- 
itational energy \Eg\ as functions of the metric function 
g2A(_R) ^']^^ surface. Compared to the corresponding 
frequencies scaled by the neutron star mass as functions 
of M/R, while the universal behavior of the wi- and w//- 
mode frequencies remain about the same, the w/2-mode 
frequencies scaled by the total gravitational energy as 
functions of e^"^^^^ become much more universal. 

Besides the above two kinds of universal scalings, it is 
possible to obtain alternative scalings. One particularly 
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To be more quantitative and make it easier to compare 
with future studies, we summarize in Table II several neu- 
tron star properties, the corresponding eigen-frequencies 
and damping times for the EOSs considered in this work. 

Finally, in Fig. 12 we show both the frequencies and 
the damping times of the w/-, w//-, w/2-, and w/3-modes 
for canonical neutron stars of mass M = IAMq. In order 
to evaluate their dependence on the symmetry energy, we 
used all four symmetry energy functionals shown in Fig.l, 
namely x = 0,-0.5,-1 and the APR, that are within the 
constraints set by the heavy-ion reaction data. It is seen 
that the various modes have different dependences on 
the symmetry energy. While the damping time of the wj 
mode increases with the increasing frequency as the sym- 
metry energy becomes softer (from a; = — ltoa; = Oas 
shown in Fig. 1), the opposite behaviors are observed for 
the W12, Wis a-iid wn modes. Thus, simultaneous studies 
of multiple w-modes will be useful in understanding the 
imprints of symmetry energy on the w-modcs of gravita- 
tional waves. It is also seen that the w/j-mode has the 
smallest frequency and damping time - the frequency of 
iti/z-mode for a neutron star of mass M = IAMq with 
the x = -1 EOS is about 1.72 kHz. 
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V. SUMMARY 

In summary, we have examined the eigen-frequencies 
of the first few axial w-modes of oscillating neutron stars 
by using an EOS with the symmetry energy partially 
constrained by the recent terrestrial nuclear laboratory 
data. Our studies indicate that the density dependence of 
the nuclear symmetry energy Egym (p) has a clear imprint 
on both the frequency and the damping time of the axial 
w-modes. We confirmed the previously found universal 
behavior of the mass-scaled eigen-frequencies. Moreover, 
we explored several alternative universal scalings of the 
eigen-frequencies. The latter scaled with the absolute 
value of the gravitational energy \Eg\ are more universal 
functions independent of the symmetry energy, especially 
for the w/2-mode. Furthermore, it is found that for the 
w/j-mode to exist neutron stars have to have a minimum 
compactness of M/R « 0.1078 independent of the EOS 
used. 



useful scaling we observed is that the real frequency u) 
of the w/j-mode scaled by the mass M varies linearly 
and independent of the EOS as a function of the metric 
function e2A(_R)^ pjg^ 21. From the lower panel of 
Fig. 3, one can see that there exists a minimum mass or 
compactness, below which the w/z-mode frequency will 
vanish, namely, there is no w/j-mode below the limit. 
Unfortunately, it is not easy to obtain this limit directly 
based on Fig. 3. However, the linear characteristic shown 
in Fig. 11 allows one to easily and accurately determine 
the minimum compactness for the existence of the wn- 
mode to be M/R « 0.1078. 
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TABLE II: Neutron star properties (central density, radius, mass, compactness and gravitational energy), and eigen-frequencies 
and damping times of the wi- and wn- modes for the EOSs considered in this paper. 



EOS 




R{km) 


M{Mq) 


M/R 


\E,\{Mq) 






Vw,j {kHz) 


r^.i, (/is) 




0.7265 


12.25 


1.00 


0.121 


0.108 


8.35 


20.8 


0.67 


11.9 




0.7970 


12.16 


1.10 


0.134 


0.133 


8.25 


22.6 


1.23 


12.2 




0.8728 


12.07 


1.20 


0.147 


0.161 


8.18 


24.5 


1.75 


12.6 




0.9586 


11.96 


1.30 


0.160 


0.194 


8.11 


26.7 


2.26 


12.9 


MDIxO 


1.0595 


11.83 


1.40 


0.175 


0.232 


8.04 


29.2 


2.81 


13.3 




1.1780 


11.68 


1.50 


0.190 


0.276 


7.98 


32.2 


3.36 


13.6 




1.3196 


11.49 


1.60 


0.205 


0.327 


7.93 


35.9 


3.95 


14.0 




1.4969 


11.25 


1.70 


0.223 


0.390 


7.89 


40.8 


4.61 


14.4 




1.7760 


10.93 


1.80 


0.243 


0.469 


7.84 


48.2 


5.36 


15.0 




0.5535 


13.94 


1.10 


0.116 


0.113 


7.30 


23.0 


0.39 


13.5 




0.6140 


13.83 


1.20 


0.128 


0.137 


7.26 


24.7 


0.81 


13.8 




0.6835 


13.70 


1.30 


0.140 


0.165 


7.22 


26.5 


1.25 


14.1 




0.7630 


13.54 


1.40 


0.153 


0.197 


7.21 


28.5 


1.72 


14.4 


MDIx-1 


0.8570 


13.35 


1.50 


0.166 


0.234 


7.21 


30.8 


2.17 


14.6 




0.9700 


13.12 


1.60 


0.180 


0.278 


7.22 


33.6 


2.69 


14.9 




1.1181 


12.84 


1.70 


0.195 


0.330 


7.24 


37.0 


3.25 


15.2 




1.3220 


12.47 


1.80 


0.213 


0.395 


7.29 


41.6 


3.90 


15.5 




1.6760 


11.90 


1.90 


0.236 


0.488 


7.39 


49.5 


4.82 


15.8 




0.7880 


11.66 


1.00 


0.127 


0.114 


8.73 


20.6 


0.99 


11.5 




0.8383 


11.62 


1.10 


0.140 


0.140 


8.60 


22.5 


1.57 


11.9 




0.8895 


11.58 


1.20 


0.153 


0.169 


8.45 


24.6 


2.08 


12.2 




0.9422 


11.55 


1.30 


0.166 


0.201 


8.31 


26.9 


2.55 


12.7 


APR 


1.0019 


11.51 


1.40 


0.180 


0.237 


8.17 


29.5 


3.03 


13.1 




1.0665 


11.46 


1.50 


0.193 


0.277 


8.03 


32.4 


3.50 


13.5 




1.1370 


11.41 


1.60 


0.207 


0.322 


7.89 


35.9 


3.97 


14.0 




1.2155 


11.32 


1.70 


0.221 


0.374 


7.75 


40.1 


4.43 


14.6 




1.3140 


11.25 


1.80 


0.236 


0.432 


7.61 


45.4 


4.91 


15.2 
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